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Abstract 

The paper takles with a procedure which allow to extend some linear, wave type models cur- 
rently used in describing phenomena appearing in atmosphere to the study of nonlinear models. 
More concretely, we present a practical way to generate the largest class of (1 + l)-dimensional 
second order partial differential equations (pdes) of a given form which could be reduced to an 
imposed ordinary wave type equation. This class generalize the ordinary differential equation 
describing the equatorial trapped waves generated in a continuously stratified ocean and will 
be obtained following the Lie symmetry and similarity reduction procedures. Moreover, some 
concrete nonlinear second order differential equations will be proposed as possible candidates 
for replacing more complicated, nonintegrable systems, as the Rossby type equation. 

Keywords: Nonlinear dynamical systems, Lie symmetries, Similarity reduction procedure, 
Rossby type symmetries. 

1 Introduction 

A rich variety of complex phenomena occuring in many physical fields, including the atmospheric 
dynamics, are described by linear differential equations which allow a simple handling of the con- 
straints which appear in the system's evolution. Although, the linearized models often do not 
adequately describe the dynamics of the processes as a whole and it is very simple to shift the 
system to a region in which the linear behavior is no longer valid. This is why, in order to capture 
the real behavior, to accurately estimate and control the complex systems in all their regimes, 
nonlinear models must be defined. In this case, the linear differential equations could appear as 
approximations to the nonlinear systems, valid under restricted conditions. 

The price to be paid when nonlinearity is taken into consideration appears in the investigation 
of the exact solutions of the attached equations. There are not standard methods of solving 
nonlinear differential equations, they are usually depend on the form of the equations and on their 
particular symmetries. Many interesting nonlinear models have been proposed over the last years 
[T] and a lot of methods have been developed in order to find solutions of equations describing 
these nonlinear phenomena. Some of the most important methods [7] are the inverse scattering 
method [8], the Darboux and Backlund transformations [9], the Hirota bilinear method |10j . the Lie 
symmetry analysis [HI H21 E3] , etc. By applying these methods, many types of specific solutions 
have been obtained. For example, solitary waves or solitons, which have no analogue for linear 
partial differential equations, are very important for the nonlinear dynamical systems. 

In this paper we shall concentrate our attention to the Lie group method. It is well-known 
that this method is a powerful and direct approach to construct many types of exact solutions 
of nonlinear differential equations, such as soliton solutions, power series solutions, fundamental 
solutions [TH [TS], and so on. The existence of the operators associated with the Lie group of 
infinitesimal transformations allows the reduction of equations to simpler ones. The similarity 
reduction method for example is an important way of transforming a (1 + 1)— dimensional pde into 
an ordinary differential one. We shall concretely consider the inverse symmetry problem [16J and 
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we shall generate the largest class of second order (1 + 1)— dimensional pdes which generalize the 
ordinary, wave type, differential equation describing the equatorial trapped waves generated in a 
continuously stratified ocean. Practically, four types of waves appears in this case and have to 
be found among the solutions of the equation: Kelvin waves, Rossby waves, inertia-gravity waves 
and mixed Rossby-gravity waves. In the Boussinesq approximation and on an equatorial /3-plane, 
the equation which describe the m-th oscillation mode of the wave's vertical velocity <j> m (z) on the 
direction z has the form |23| : 

d 2 cp m {z) Nj(z) 

dz 2 + "^-iW = C 1 ) 

One consider for the velocity the boundary conditions: 

0m( z = ~~ H) = (ocean floor) 
§mk z = 0) = (ocean surface) 

In the equation ([T|), C m is a constant and N m (z) represents the "buoyancy" frequency. Mea- 
surements made during El Nino events [23] show that the buoyancy frequency N m (z) has strong 
variations with the water depth close to the surface and practically vanishes for higher depths. We 
notice that in the first case (at the surface, z G [0,300] m) one can aproximate N m (z) with an 
averaged value around N(z) = 2 • 10 _4 m • s _1 So, the equation ([1]) can be linearized in one of the 
following forms: 

<j)(z) = 0;z> 300 (2) 

<j){z) + k 2 <\>{z) = 0; z e [0, 300]; k = — = const. (3) 

In this paper, we shall consider the two previous wave type equations and we shall see how 
they can be extended towards (1 + 1)— second order differential equations with the same group of 
symmetries as the initial ordinary wave type equations have. 

The outline of this paper is as follows: after this introductory notes, in Section 2, we shall obtain 
the general determining system for a chosen class of (1 + 1)— dimensional models. The system will be 
generated by using the Lie symmetry approach and by asking for an imposed form of the similarity 
reduction equation. More exactly, we shall generate a class of (1 + 1)— second order differential 
equations which by similarity reduction come to the wave forms ([2]) and, respectively, Q. The 
general results of the second section will be particularized in Section 3, when concrete examples of 
two dimensional equations with similar solutions as the ordinary wave equations ([2]) and will be 
generated. Moreover, we shall compute the form of the second order partial differential equation 
which admit an imposed form of symmetry, specific for the two dimensional Rossby type equation. 
So we shall be able to replace the study of this last strongly nonintegrable equation with a simpler 
class of equations observing similar symmetries. Some concluding remarks will end the paper. 



2 Determining equations for the Lie symmetry group 

Let us consider the class of general dynamical systems described in a (1 + l)-dimensional space 
(x,t) by a second order partial differential equation of the form: 

u t = A(x, t)u2x + B(x, t)u x + C(x, t)u 44> 0(x, t, u, u x ,ut,U2 X ) = (4) 

Our aim is to select from the general dynamical systems described by @) the class of differential 
equations which admit a similarity reduction to wave type equations of the form (J2|) and Q. As 
feed-back, the solution of the wave equations will be used in order to obtain a solution for @. The 
procedure that will be followed firstly implies to obtain the system of determining equations for 
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the Lie symmetry group of Then, an additional system of partial differential equations will be 
generated by imposing that possess a reduced similarity equation of the wave type. Finally, this 
last system and the Lie determining equations will be solved and the coefficient functions A(x,t), 
B(x,t), C(x,t) will be obtained. 

In this section we shall apply the Lie symmetry approach for the equation Let us consider 
a one-parameter Lie group of infinitesimal transformations: 

x = x + e£(t, x, u), t = t + ecp(t,x,u), u = u + erj(t,x,u) (5) 

with a small parameter e <C 1. The Lie symmetry operator associated with the above group of 
transformations can be written as: 

d d d 

U(x, t, u) = ip{x, t,u)—+ t,u)— + 7](x, t, u)— (6) 

The second order equation Q(x,t,u, u x ,ut,U2 X ) = of the form ^) is invariant under the action of 
the operator (0) if and only if the following condition [TT] is verified: 

£/ (2) (^) b=o=0 (7) 

where is the second extension of the generator ([6]). A concrete computation shows that the 
coefficient functions from (jU) and ([6]), A(x, t), B(x, t), C(x, t), ip(x, t, u), £(x, t, u), rj(x, t, u), must 
satisfy the equation: 

(<pA t + iA x )u 2x + (tpBt + £B x )u x + LpC t u + £C x u + Crj + Brf - rf + Ar] 2x = (8) 

The coefficient functions T) x , rf, r] 2x appear in the process of extension of U towards and their 
general expressions are given in . Using these expressions in ([$]) and asking for the vanishing of 
the coefficients of each monomial in the derivatives of u(t,x), we obtain the following differential 
system: 

¥>* = 0; v u = 0; = 0; V 2u = 0; (pA t + £A X + Atp t - 2A£ X = 0; 

- ipB t - iB x + B£ x - & - Bcp t - 2Ar ]xu + A£ 2x = (9) 

-ipC t u - £C x u -Cr]- Br] x + i] t + Crj u u - (p t Cu - Arj 2x = 

The first four equations of the system ([9|) lead, for coefficient functions ip(x, t,u),£ (x, t, u),rj(x, t, u), 
to the following reduced dependences: 

if = <p(t), £ = C(X, t), 7] = M(x, t)u (10) 

Consequently, the remaining equations of (|9j) become: 

V A t + iA x + Asp t - 2Ai x = 0; 

- ipB t - iB x + B£ x - £ t - B<p t - 2AM X + A£ 2x = (11) 



-ipCt - £C X - BM X +M t - <p t C - AM 2x = 

with 6 unknown functions: A(x,t), B(x,t),C(x,t) provided by the evolutionary equation and 
ip(t), £(x,t), r}(x,t,u) introduced by the symmetry group of transformations ([5]) and described by 
the relations (jlOp . 
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3 Similarity reduction procedure 

Let us consider now the similarity reduction procedure. In this section, some particular choices for 
the system (|lip will be considered. We shall find equations describing concrete dynamical systems 
which admit reduction through the similarity procedure to ordinary wave type equations of the form 
([2]) and ©. For the moment, we restrict the forms (fTUj) of the infinitesimals ip(t),£(x,t),r](x,t,u) 
to the following separable expressions: 

ip = <p{t), £ = £(x,t) = £i(s)£ 2 (t), rj = M{x,t)u = M 1 {x)M 2 (t)u (12) 

The Lie operator @ becomes: 

U(x,t,u) = ip(t)-+Z 1 (x)Z 2 (t)— + M 1 (x)M 2 (t)u— (13) 

The general expressions of the invariants could be found if we should consider the characteristic 
equations associated with the new generator (1131) . These equations are: 



dt dx du 

Ri) = £i(aO£ 2 (*) = M 1 {x)M 2 (t)u ( } 

By integrating the previous equations, two invariants are obtained with the following expressions: 

In the similarity reduction procedure two similarity variables have to be considered: 

h = z, h = <f>(z) (16) 

The invariants (fT6j) allow us, by an appropriate change of coordinates {u, x,t} — > {z,(p(z)}, to 
reduce the initial (1 + 1) dimensional equation to an ordinary differential equation of the form: 

n'[z,<i>(z),j)(z),...}=0 (17) 
3.1 Homogeneous wave type equation 

Our aim is now to select from the general dynamical systems described by (dJ the class of differential 
equations for which the equation (|17p can be reduced at a wave type equation of the form ([2]): 

1 = o (f)(z) = az + b (18) 

where a and b are arbitrary constants. 

The previous solution, written in terms of the initial variable (x,t), leads to the following form 
of the solution u(x,t) of (jlj): 



u(x, t) 



aexp ( / dx — f ^ 2 dt ) + b 

Cx(x) J <p(t) 



(MM ™ 

p v e 2 (*) J Mx) 



(19) 



For convenience reasons, we shall impose the following relations to be valid: 
Z 2 (t) = M 2 (t) 

with q, v arbitrary constants. 
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In terms of notations (|20p . the infinitesimals (|12p and the solution (|19p become: 



u(x, t) 



'Pit), £ 



^— 1i 



(21) 
(22) 



'P(x)' ' ' P(x) 
[aexp (P(x) — qt) + 6] exp (vR(x)) 

The solution (j22[) must verify the equation ([3]) which describes the analyzed model. This condition 
generates a differential system of the form: 

= q + 2vA(x,t)R(x)P(x) + v 2 A(x , t)[R(x)] 2 + A(x,t)P(x) + A(x,t)[P(x)] 2 + 

+vA(x, t)R(x) + 5(x, t)P(x) + vB(x, t)R(x) + C(x.t) 
= u 2 ^(a;,t)[fl(x)] 2 + vA(x,t)R(x) + vB(x,t)R(x) + C(x,t) (23) 

For an unitary analysis, it is necessary to describe the differential system obtained in 

the previous subsection, in terms of the functions P(x) and R(x) introduced by (p0|) . Using the 
expressions (|2ip we obtain the following differential system: 



ipA t P 2 + qipA x P x + ip t AP 2 + 2qipAP 2x 







pB*i* + ^P.Pj + q<pBP 2x P 2 + + p t flj£+ 

+ 2vq<pAR 2x P* - 2vqipAR x P 2x P 2 + q^AP 3x P 2 



(24) 



2g^P x P4 = 

</?C t P 4 + gpC x i£ + qv^BR 2x Pl - qvipBR x P 2x P 2 + (^CP 4 + vq V AR 3x P^- 

+vq<pAR x P 3x P 2 - 2vqy?AR 2x P 2x P 2 + +2^,4^^ - qv^ t R x P^ = 

Conclusion: Our problem is to find the class of (1 + 1) evolutionary equations of type (HD 
which could be reduced by the similarity approach to an ordinary wave type equation. Solving this 
problem is equivalent with searching the solutions of the system described by equations (|23|) and 

HMD- 
Remark 1: The system (j23p ~(|24 |) can be solved following two paths: (i) by choosing a concrete 
dynamical system, that is to say concrete expressions for the functions A(x,t), B(x,t), C(x,t) and 
trying to find out if this equation admits or not solution of the type f)22|) . Now the unknown 
functions of the system are <p(t), P(x), R(x) defined by (p0|) : (ii) by considering A(x,t), B(x,t), 
C(x,t) as unknown functions and by choosing <p(t), P(x), R(x). This is the way we shall follow in 
the next section. 

Remark 2: In the case (ii) the general solutions obtained by computational way can be 
expressed as: 

'qt - P(x)' 



A(x,t) 
B(x,t) 
C(x,t) 



-2F 



qt-P{x) 

q 



exp [Cr(ac)] 

\[P(x)} 2 + vR(x)P(x) + lP(x)) exp [-G(x) 



qt-P{x) 

q 



P(x) 

R(x)P(x) + P(x)(-R(x) + R(x)(vR(x) + P(x)))] exp^ G ^ +qR(x) 



where 



G(x) 



P(x) 



2(D^)(P)(a)ip ( P(a)+<lt-P(x) \ g + [jD ( P )( a )]2 D{lp) ( P(a)+ q t-P(x) 



D(P)(a)cp 



P{a)+qt-P(x) 

q 



da 



(25) 



(26) 



These solutions are valid for arbitrary constants q, v and for an arbitrary function F 



qt-P(x) 

q 
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3.2 Harmonic oscillators 



Let us consider now that the similarity reduction equation is an ordinary oscilator type equation 
of the form ([3|): 

^ + *V(*)=0 (27) 

It has the solution: 

4>(z) = asin(kz) + bcos(kz) (28) 

Here k, a and b are arbitrary constants. 

To writte down the previous solution in terms of the initial variable u(x, t) means that the 
solution of (jH) should have the form: 

u(x,t) = aexp ( f Ml ^ X ) dx] sin ( Vfcexp ( [ -^--dx - [ QQ-dtX) + (29) 

\w) J £iw / V \J J // 

+6exp ( —tt-t- [ - r- | cos ( \/&exp ( / dx — f } dt 
V V fi(a?) / V U CiW J tp(t) , 



For convenience reasons, we shall impose again the following relations to be valid: 

e,w _ . W) _ „. t i dl s PW , / ^M,, 3 (30) 



*>(<) " &(«) J fiW " 1 '7 

with g, u arbitrary constants. 

In terms of notations (|3(jp . the infinitesimals (|12p and the solution (|29|) become: 

^ = ¥>(t), £ = ?t— . V = QV . u (31) 
P[x) F{x) 



u(x,t) = [asin(/cexp (P(x) — qt)) + 6cos(/cexp (-P(x) — <#))] exp (yR(x)) 

(32) 

The solution (j32|) must verify the equation (J3J) which describes the analyzed model. This condition 
generates the vanishing of the coefficient function A(x,t) and two other differential equations of 
the form: 

A(x,t) = 

q + B(x,t)P(x) = (33) 
vB(x,t)R(x) + C(x,t) = 

For an unitary analysis, it is again necessary to describe the general differential system 
obtained in the previous section, in terms of the functions P(x) and R(x) introduced by ([30]) . 
Taking into account the equations (|3ip . we obtain the following differential system: 

= vB t Px+qvB x Pl + q<pBP 2x Pl+q<p t Pl + VtBPt 
= i P C t P x i + qipC x P x i + qvLpBR 2x P x i -qvipBR x P 2x P x 2 + 

+^ t CPt - qv^tRxPl (34) 

The system (|33|) - (|34|) can be solved following two paths: {%) by choosing a concrete dynamical 
system, that is to say concrete expressions for the functions B(x,t), C(x,t) and trying to find out 
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if this equation admits or not solution of the type ()32j) . Now the unknown functions of the system 
are <p(t), P(x), R(x) defined by ([30]) ; (ii) by considering B(x,t), C(x,t) as unknown functions and 
by choosing tp(t), P(x), R(x). 

This second case is the way we are interested in to follow and, in this case, the general solutions 
obtained by computational way can be expressed as: 

B (,, ( )--Z« (x, ( )-!**5> (35) 
P(x) P(x) 

or in terms of the coefficient functions <p(t), £(x,t), M(x,t) which appear in the general Lie sym- 
metry operator (??), in the equivalent forms: 

B(x,t) = ^H, C M - MM 



3.3 Rossby type symmetries 

The equation for coupled gravity, inertial and Rossby waves in a rotating, stratified atmosphere 
using the /3-plane approximation (which simplifies the spherical geometry whilst retaining the 
essential dynamics) and the Boussinesq approximation which filters out higher frequency acoustic 
waves can be written in (2 + 1)— dimensions in the form 



d ( d 2 &_\ _ _ R duy 

dt \dx 2 + dy 2 ) U y-~P dx (36) 

As we mentioned, this equation describes the coupling between the inertial, the gravity and the 
Rossby waves, but also the shallow water in an ocean of depth H. It was proven [25] that reducing 
the model to (1 + 1)— dimensions, (x,t), it admits some very simple Lie symmetries of the form: 

<p(t) = cb + ci, £(x,t) = cx + f(t), r](u) = —3cu (37) 

with f(t) arbitrary function and cq, c arbitrary constants. Despite this fact, it is still difficult to 
find explicite solutions for the equation (|36p . This is why, we shall consider another approach: we 
shall impose the Lie symmetries (|37]1 to our general equation @ and we shall try to find the class of 
equations which observe them. This means that we have in fact to impose the Rossby symmetries 
(|37p to the system (|lip . It will take the form: 



(ct + d)A t + {cx + c 2 )A x + 3cA = 
(ct + cx)B t - (cx + c 2 )A x -2cB = 
-(ct + d)C t - (cx + c 2 )C x - cC = 



with the unknown functions A(x,t), B(x,t), C(x,t). 
This system admits the solutions: 



A(x,t) 
B(x,t) 
C(x,t) 



F(x(ct + ci) - 


P2*) 


F (x<p(t) - 


c 2 t) 


(ct + Cl) 3 








G (x(ct + ci) - 


c 2 t) 


G (xtp(t) - 


c 2 t) 


(ct + ci) 2 




W{t)] 2 




H(x(ct + ci) - 


c 2 t) 


H (xip(t) - 


c 2 t) 


ct + Cl 




<p(t) 





with F, G, H arbitrary functions of their arguments. These expressions give us equations of the 
form Q which are equivalent from the point of view of their symmetries with the Rossby equation. 
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4 Conclusions 



The problem of finding exact solutions for nonlinear differential equations plays an important role in 
the study of nonlinear dynamics. There are many ways of tackling with it. One of them is based on 
the Lie symmetry method. This method supposes to find the symmetries of the system and, on this 
basis, to try to determine the general or some particular solutions of the equations. There is a direct 
approach in which the symmetries of a given equation are obtained, but also an inverse problem has 
been formulated [16J. A step forward for this latter approach is represented by the use of similarity 
reduction, a procedure which allows the reduction of the number of degrees of freedom and, by 
that, simplifies the problem of solving the equation. This paper used this approach and determined 
a class of (1 + 1) dimensional second order differential equations which can be reduced to ordinary 
wave-type equations with simple solutions. Using the Lie symmetry and the similarity reduction 
procedures, some particular cases of the equation arise as good candidates of equations which 
could be used as generalization of the linear wave type equations describing complex atmospheric 
phenomena. Moreover, following our method, we were able to write down the solutions of these 
equations, solutions which otherwise could be derived by computational methods, but in a very 
complicated form. Another interesting results of our paper consisted in the fact that a complicated, 
nonintegrable equation, the Rossby equation, could be replaced by another, simpler equation, which 
have similar symmetries. The paper is important both by these results, but also as a methodological 
approach in finding exact solutions through similarity reduction procedure. We have shown how, 
starting from a particular form of solution for the reduced equation we could recover the solution 
of a most complicated problem, defined in a space with more than one dimensions. We tackled 
out a particular case, looking only for linear solutions of the reduced equation and considering that 
the coefficient functions appearing in the symmetry operators are separable. The problem can be 
extended for other cases, too. 

Acknowledgements: This work was supported by CNCSIS -UEFISCSU, project number PNII 
- IDEI code 418/2008. 
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